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Abstract. We have analyzed the projected galaxy distri- 
butions in a subset of the ENACS cluster sample, viz. in 
those 77 clusters that have z < 0.1 and Raco > 1 and 
for which ENACS and COSMOS data are available. For 
20 % of these, the distribution of galaxies in the COSMOS 
catalogue does not allow a reliable centre position to be 
determined. For the other 62 clusters, we first determined 
the centre and elongation of the galaxy distribution. Sub- 
sequently, we made Maximum-Likelihood fits to the distri- 
bution of COSMOS galaxies for 4 theoretical profiles, two 
with 'cores' (generalized King- and Hubble-profiles) and 
two with 'cusps' (generalized Navarro, Frenk and White, 
or NFW, and de Vaucouleurs profiles). 

We obtain average core radii (or characteristic radii for 
the profiles without core) of 128, 189, 292 and 1582 kpc 
for fits with King, Hubble, NFW and de Vaucouleurs pro- 
files respectively, with dispersions around these average 
values of 88, 116, 191 and 771 kpc. The surface density of 
background galaxies is about 4 10~ 5 gals arcsec -2 (with a 
spread of about 2 10~ 5 ), and there is very good agreement 
between the values found for the 4 profiles. There is also 
very good agreement on the outer logarithmic slope of 
the projected galaxy distribution, which is that for the 
non-generalized King- and Hubble-profile (i.e. fixing = 
P Hubble — Ij with the corresponding values for the two 
other model-profiles). 

We use the Likelihood ratio to investigate whether the 
observations are significantly better described by profiles 
with cusps or by profiles with cores. Taking the King and 
NFW profiles as 'model' of either class, we find that about 
75 % of the clusters arc better fit by the King profile than 
by the NFW profile. However, for the individual clusters 
the preference for the King profile is rarely significant at 



a confidence level of more than 90 %. When we limit our- 
selves to the central regions it appears that the signifance 
increases drastically, with 65 % of the clusters showing a 
strong preference for a King over an NFW profile. At the 
same time, about 10 % of the clusters are clearly better 
fitted by an NFW profile than by a King profile in their 
centres. 

We constructed composite clusters from the COSMOS 
and ENACS data, taking special care to avoid the creation 
of artificial cusps (due to ellipticity) , and the destruction 
of real cusps (due to non-perfect centering). When adding 
the galaxy distributions to produce a composite cluster, 
we either applied no scaling of the projected distances, 
scaling with the core radii of the individual clusters or 
scaling with r2oo, which is designed to take differences 
in mass into account. In all three cases we find that the 
King profile is clearly preferred (at more than 95 % confi- 
dence) over the NFW profile (over the entire aperture of 5 
core-radii). However, this 'preference' is not shared by the 
brightest (M b . -18.4) galaxies. We conclude that the 
brighter galaxies are represented almost equally well by 
King and NFW profiles, but that the distribution of the 
fainter galaxies clearly shows a core rather than a cusp. 

Finally, we compared the outer slope of the galaxy 
distributions in our clusters with results for model cal- 
culations for various choices of fluctuation spectrum and 
cosmological parameters. We conclude that the observed 
profile slope indicates a low value for fio- This is consis- 



of the individual clusters. 
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1. Introduction 

Until fairly recently, the projected galaxy density in rich 
galaxy clusters was generally described by King or Hub- 
ble profiles. In these profiles, the logarithmic slope of the 
mass distribution is essentially zero near the cluster cen- 
tre. The core radius which is the characteristic scale of 
the distribution, was sometimes also regarded as the dis- 
tance which more or less separates dynamically distinct 
regions in a cluster. From the kinematics of the galaxy 
population it appears that in clusters the relaxation time 
is significantly shorter than the Hubble time only in the 
very central region within at most a few core radii (see 
e.g. den Hartog and Katgert 1996). 

The concept of cores in clusters has been seriously 
challenged, on observational grounds (e.g. Beers & Tonry 
1986) and as a result of numerical simulations. Navarro, 
Frenk and White (1995, 1996) found e.g. that the equilib- 
rium density profiles of dark matter halos in universes with 
dominant hierarchical clustering all have the same shape, 
which is essentially independent of the mass of the halo, 
the spectrum of initial density fluctuations, or the values 
of the cosmological parameters. This 'universal' density 
profile (NFW profile hereafter) does not have a core, but 
has a logarithmic slope of -1 near the centre which, at 
large radii, steepens to -3, and thus closely resembles the 
Hcrnquist (1990) profile except for the steeper slope of the 
latter at large radii of -4. 

Navarro, Frenk and White (1997) argue that the ap- 
parent variations in profile shape, as reported before, can 
be understood as being due to differences in the charac- 
teristic density (or mass) of the halo, which sets the linear 
scale at which the transition of the flat central slope to 
the steep outer slope occurs. They also argued that the 
existence of giant arcs in clusters requires that the mass 
distributions in clusters does not exhibit a flat core in the 
centre. In other words: if clusters have cores, the lensing 
results require that the core radii are very small, at least 
quite a bit smaller than the values usually quoted. 

It is not clear that galaxy clusters should have cores; 
after all, the dynamical structure of galaxy clusters is quite 
different from that of globular clusters, for which Michie 
& Bodcnhcimer (1963) and King first proposed density 
profiles with cores, in particular the King profile (see e.g. 
King 1962). On the other hand, the X-ray data for clus- 
ters are quite consistent with the existence of a core in 
the density distribution. More specifically, it was argued 
recently by Hughes (1997) that the NFW profile would 
induce a temperature gradient. The existence of such a 
gradient in the Coma cluster can be excluded at the 99% 
confidence level. Similarly, the galaxy surface density in 
clusters is generally found to be consistent with a King 
profile. For galaxy clusters, little use has been made of 
the de Vaucouleurs profile to describe the galaxy density, 
even though the latter was found to arise quite naturally 



in N-body simulations of the collapse of isolated galaxy 
systems (e.g. van Albada 1982). 

In view of the claimed universality of the NFW pro- 
file found in the simulations, it seems useful to have a 
closer look at the projected distribution of the galaxies 
in clusters. After all, the NFW profile refers to the total 
gravitating mass, and it is not obvious that the galaxy dis- 
tribution should have exactly the same shape as the distri- 
bution of total mass; although in numerical experiments 
no strong biasing between dark and luminous matter in 
clusters was seen (e.g. van Kampen 1995). In this respect, 
it is noteworthy that Carlberg et al. (1997) find that the 
combined galaxy density profile of 16 high- luminosity X- 
ray clusters at a redshift of «0.3 closely follows the NFW 
profile. More precisely, the logarithmic slope in the central 
region is consistent with the value of -1 of the NFW and 
Hcrnquist profiles, while the outer slope is consistent with 
both -3 (the NFW value) and -4 (the Hcrnquist value) . 

The outer slope of the density profile was found by sev- 
eral authors (e.g. Crone et al. 1994, Jing et al. 1995, and 
Walter and Klypin 1996) to reflect the details of the for- 
mation scenario, and in particular the value of the density 
parameter of the universe. In addition, this slope is un- 
likely to be constant in time but is expected to get steeper 
with decreasing redshift. For that reason, it is important 
to study both the characteristics of the density profiles of 
rich clusters and their dependence on redshift. 

In this paper, we investigate the projected galaxy dis- 
tributions for a sample of 62 rich and nearby (z 0.1) 
clusters. These clusters are taken from the volume-limited 
ENACS (ESO Nearby Abell Cluster Survey) sample of 
R-aco > 1 clusters (see e.g. Katgert et al. 1996 (paper I), 
Mazurc ct al. 1996 (paper II), Biviano et al. 1997 (paper 
III), Adami ct al. 1998 (paper IV), Katgert et al. 1998 (pa- 
per V) and de Theije & Katgert 1998 (paper VI)). In § 2, 
we first describe the sample of clusters that we used, and 
the data on which we based our analysis. In § 3, we discuss 
the results of Maximum-Likelihood fitting of profiles with 
and without a core, to the individual galaxy distributions 
taken from the COSMOS catalogue. In § 4 we discuss the 
galaxy density distribution for composite clusters (COS- 
MOS and ENACS), in which the individual clusters are 
combined. In § 5 we discuss the constraints provided by 
the outer slope of the density distributions for the param- 
eters of the formation scenario and in § 6 we present the 
conclusions. 

2. The data 

2.1. The galaxy catalogues 

In this paper we use both COSMOS and ENACS data. 
The COSMOS data, i.e. photometric galaxy catalogues 
that were obtained from automatic scanning of UK 
Schmidt Illa-J survey plates with the Edinburgh plate- 
scanning machine, were kindly provided by Dr. H.T. 
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MacGillivray. As described in paper V, the COSMOS data 
that we used arc of two kinds: the well-calibrated part 
around the Southern Galactic Pole (the so-called EDSGC, 
see e.g. Heydon-Dumbleton et al. 1989), and the slightly 
less well-calibrated remainder. In paper V, we compared 
the COSMOS and ENACS photometry and found only 
a small difference in the calibration quality of the two 
COSMOS subsets. We did not find evidence for system- 
atic magnitude offsets between the two COSMOS subsets, 
nor for differences in completeness. From a comparison 
with the ENACS galaxy catalogues, which are based on 
completely independent scanning with the Leiden Obser- 
vatory plate-measuring machine, we concluded that the 
COSMOS catalogue is 90% complete to a nominal limit 
of m^. « 19.5. 

The galaxy catalogues that resulted from the ENACS 
spectroscopic survey are described in papers I and V. Very 
briefly, redshifts were obtained for 5634 galaxies in 107 
ACO cluster candidates, mostly in the southern hemi- 
sphere. As shown in paper II, the ENACS allowed us to 
construct a complete, volume-limited sample of 128 rich 
(R-aco > 1) clusters, out to a redshift of 0.1, when we 
combine ENACS data for about 80 clusters with literature 
data for about 50 clusters. In paper V, we have shown the 
magnitude distributions of the galaxy samples for which 
we obtained ENACS spectroscopy. Comparison with the 
magnitude distributions of the COSMOS galaxies shows 
the ENACS samples to be more or less complete to 
between about 18 and 19, with quite a few redshifts for 
fainter galaxies, viz. down to m^. « 19.5. 

2.2. The cluster sample 

The overlap of the ENACS dataset with that part of the 
COSMOS dataset that was available to us yields a sample 
of 77 clusters. These clusters were not selected according 
to particular criteria, and we therefore expect them to be 
a representative subset of the total ENACS sample. To 
this sample, we have added the cluster A2721, for which 
wc use redshift data from Collcss & Hcwett (1987), Colless 
(1989) and Teague et al. (1990). 

Because we want to study density profiles, we are re- 
luctant to use clusters with clear signs of substructure, as 
for the latter a scale length and the inner and outer slopes 
of the density profile do not have well-defined meanings. 
In addition, it is not at all trivial to define a meaningful 
centre for systems that appear irregular in projection. 

In order to give a visual overview of the galaxy distri- 
butions in the 77 clusters in the sample, we have used the 
COSMOS data to produce adaptive-kernel maps for all 77 
clusters; these are shown in Figs. 1 and 2. 

The contour maps in these figures are based on sur- 
face densities calculated as the sum of normalized gaus- 
sian kernel functions at the positions of the galaxies with 
dispersions Oi that are adapted to the local galaxy density 
(e.g. Silverman 1986). In Figs. 1 and 2 , it is fairly easy 



to distinguish clusters with single density peaks, clusters 
with multiple density peaks that can be clearly separated 
and/or identified with systems at different redshifts, and 
irregular clusters. By using different symbols for systems 
with different redshifts (see papers I and IV), it is possible 
to recognize clusters which consist of several 'single' peaks 
at different redshifts. About two-thirds of the clusters in 
our sample show a single peak, while about 10% was con- 
sidered irregular. We have defined the studied area for 
each cluster as the maximal area without a clear sign of 
substructure and with a single density peak. Afterwards, 
we have checked that these areas comprise more than 5 
King core radii (i.e. with a radius greater than 500 kpc). 

In Tab. 1, we give our verdict on the character of the 
galaxy distribution, as well as some other information, for 
the 62 clusters in Fig. 1 for which a central position could 
be determined reliably (see § 3). The 15 clusters for which 
no reliable position could be obtained are shown in Fig. 2; 
those could not be used in the following analysis. It should 
be noted that Tab. 1 has 2 entries for A0151, while the 
cluster A3559 (in the region of the Shapley concentration) 
was discarded because it did not show a clear maximum 
in the galaxy distribution. Note also that our verdict 'reg- 
ular' in Tab. 1 does not imply that the cluster does not 
have substructure, as projection may hide substructure 
along the line-of-sight. 

The mean redshift of the systems in Tab. 1 is about 
0.075 (practically equal to that of the total ENACS sam- 
ple). In the following analysis we will adopt Ho = 100 
km/s/Mpc and qq = 0. The parameters will be determined 
with robust statistical estimators (Beers et al. 1990). 

3. The density profiles of the individual clusters 

3.1. The Maximum-Likelihood fits 

We derived the characteristics of the projected galaxy 
distributions for the individual clusters using Maximum- 
Likelihood fits (hereafter MLM fits, see e.g. Sarazin 1980 
for a description of the method). We have made MLM fits 
for different types of profile as follows. Define a(x,y) as 
the theoretical, normalized density profile with which one 
wants to compare the data. Note that we scaled the am- 
plitude of the model to reproduce the observed number of 
galaxies. 

The probability that this assumed profile 'produces' a 
galaxy in position (xfc,yfc) is thus a(x k , y k )- Consequently, 
the combined probability that the assumed profile will 
'produce' galaxies in the positions (~x k ,y k ) (with k=l...N) 
that they actually have is: 

N 

L= [] o-(x k ,y k ) 
fe=i 

The model-parameters which produce the best fit be- 
tween the data and the model are found from a maxi- 
mization of the likelihood L. The combined probability, 
L, for all galaxies to have their actual positions assum- 
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Fig. 1. Adaptive-kernel maps of projected galaxy density from the COSMOS catalogue. Coordinates are in arcsec with respect 
to the cluster centre. Different groups of redshift are indicated as follows: nearest group on the line of sight: circles, 2nd group: 
squares, 3rd: triangles, 4th: crosses, 5th: asterisks and 6th: stars. We have indicated the area we use for the 29 clusters involved 
in the COSMOS composite clusters. For A1069, we use the whole area. 
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Fig. 1. continued 
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Fig. 1. continued 



ing the model to be correct, is always less than unity. For 
that reason, in practice, one usually minimizes the positive 
parameter 21n(L) w.r.t. the parameters of the assumed 
profile. This yields the values of the parameters for which 
the model is most likely to have generated the N galaxies 
at their observed positions. The great advantage of this 
method is that it does not require the data to be binned. 
In addition, it uses all information that is available. 

The minimization of 21n(L) was done with the two 
minimization methods in the MINUIT package: SIM- 
PLEX (Nelder & Mead, 1965) and MIGRAD (Fletcher 



1970). We use the same stategy as described in paper IV 
to obtain convergence. Viz., we use SIMPLEX to approach 
the optimal values and MIGRAD to refine those and get 
reliable error estimates. When MIGRAD did not converge, 
we adopted the SIMPLEX value without error estimates. 

We used the following model profiles for the 2D galaxy 
distributions: 

- generalized King model: <j(r) = co( i + tj_y2 ) + &b, 

- generalized Hubble model: a(r) = er ( ppz: ) 213 + a b, 

- generalized NFW model: er(r) = g ( j_(i\j_yi Y + 
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Fig. 2. Adaptive-kernel maps of projected galaxy density from the COSMOS catalogue for the 15 clusters for which the centre 
could not be determined reliably. Coordinates are in arcsec with respect to the cluster centre. The different groups of redshift 
are characterized by a different symbol: nearest group on the line of sight: circles, 2nd group: squares, 3rd: triangles, 4th: crosses, 
5th: asterisks and 6th: stars. 



Note that the 2-D model profile that we refer to as 
NFW is not an exact 2-D version of the 3-D profile de- 
scribed by NFW. Instead, it is a pseudo NFW profile, with 
a relation between 3-D and 2-D outer logartithmic slopes 
as for the King and Hubble profiles, but with an inner 
logarithmic slope equal to that of the 3-D NFW profile. 

at r \0.25 

- generalized dc Vaucoulcurs model: a(r) = a e p( r c > + 

Ob- 



As the actual clusters are not necessarily axi-symmetric, 
the models have 7 free parameters: two parameters for the 
position of the centre (x , yo), two parameters to describe 
deviations from symmetry (ellipticity e and position an- 
gle </>), two parameters that specify the profile (r c and (3) 
and, finally, the background density cr b (assumed constant 
within the aperture of each cluster) . 

In principle, an MLM solution could have been made 
for all 7 parameters simultaneously. However, we have sep- 
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Table 1. Parameters of the 62 clusters for which a reliable cen- 
tre position was obtained. Col.(l) cluster name, col. (2) ENACS 
redshift of the group, col (3) type of galaxy distribution (0: reg- 
ular single peak, 1: multimodal, 2: irregular), cols. (4) and (5) 
centre position, col (6) ellipticity (when available) and col (7) 
direction of the major axis (if ellipticity available and not equal 
to 0), counted anti-clockwise from 8=0. 



ACO z type a 5 c (p 

2000.0 



16 



25 



45 
70 
4.3 
25 





-2 



28 
52 
20 



27 



5 
35 
08 



30 
7.3 



50 



49 
61 
41 



35 
35 



28 
15 



arated the solution for the central position and the elon- 
gation and position angle, from the solution of the back- 
ground density and the profile parameters r c and /3. More 
precisely: we used a King profile with r c = 100 kpc, (5 — 
1.0, and we assumed a value for Ob of 3 10 -5 galaxies per 
square arcsec to make an MLM fit for x , yo, e, and 4>. As 
mentioned previously, this fit did not converge for 15 clus- 
ters, viz. for A0295, A0303, A0420, A0543, A2354, A2361, 
A2401, A2569, A2915, A2933, A3108, A3264, A3354, 
A3559 and A3921. These clusters could therefore not be 
used for MLM fits for the four profiles. 

The average uncertainty in the central position as 
given by the fits is about 30 kpc. This is not too differ- 
ent from, but somewhat smaller than the measured offsets 
between fitted positions and literature X-ray centres or lit- 
erature positions of cD galaxies as given in Tab. 2. Except 
for the clusters A0514, A3009 and A3128, which are ei- 
ther double-peaked or irregular, and apart from A3825 all 
offsets are smaller than 100 kpc, and the biweight average 
measured offset between fitted positions and X-ray or cD 
positions is 51 kpc. 



Table 2. Literature positions for 17 clusters with either X-ray 
or cD centre position. Col.(l) cluster name, cols. (2) and (3) 
literature centre, col. (4) distance between literature and fitted 
centre in kpc, col. (5) type of literature centre (X or cD) 
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In Tab. 3 and Fig. 3, we show a comparison with the 
more homogeneous ROSAT sample of X-ray centres (pre- 
liminary centres kindly provided by H. Bohringer). The 
biweight average offset between fitted positions and X-ray 
centres is 78 kpc. If we ignore the clusters with an atyp- 
ical offset of more than 150 kpc, we obtain 69 kpc. The 
clusters that we considered atypical are A0380 which is 
a clear double-peaked cluster, A3809 which is irregular, 
A2871 which exhibits clearly a very atypical difference of 
449 kpc, A0524 and A2799. We have checked the internal 
accuracy of the ROSAT centres by using the X-ray map 
of AO 119 (kindly provided by D. Neumann). We estimate 
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that this error is about 20" (11 kpc) for A0119. However, 
we deal we a very regular cluster and the error is probably 
underestimated compared to the other clusters. Another 
way to estimate the global error of the ROSAT centres is 
to compare with other literature positions (X-ray or cD). 
We find a mean offset of 31 kpc (see Tab. 3), greater than 
the 11 kpc obtained for A0119. 



Fig. 3. The histogram of the offsets between the X-Ray 
ROSAT and the fitted centres. The two clearly atypic values 
for A0380 and A2871 are not shown. 



Table 3. ROSAT X-ray centre positions for 28 clusters. Col.(l) 
cluster name, cols. (2) and (3) ROSAT centre, col. (4) distance 
between ROSAT and fitted centre in kpc, col.(5) distance be- 
tween ROSAT and previous literature centre in kpc (when 
available) . 



ACO 


a 


<5 


dist. 


dist. ROSAT/Lit 




2000.0 


(kpc) 


(kpc) 


0013 


00:13:38.23 


-19:30:08 


75 




0119 


00:56:14.60 


-01:15:44 


47 


40 


0151 


01:08:50.47 


-15:24:33 


18 


4 


0367 


02:36:39.34 


-19:23:09 


94 




0380 


02:44:05.69 


-26:11:17 


501 




0524 


04:57:57.19 


-19:43:58 


150 




0978 


10:20:28.72 


-06:31:14 


75 




1069 


10:39:44.81 


-08:41:01 


42 


22 


2426 


22:14:32.38 


-10:22:12 


61 


25 


2717 


00:03:11.14 


-35:55:21 


74 




2734 


00:11:20.28 


-28:51:31 


44 


21 


2755 


00:17:34.90 


-35:11:09 


97 




2764 


00:20:34.10 


-49:13:40 


72 




2799 


00:37:29.47 


-39:06:18 


152 




2871 


01:07:49.37 


-36:43:44 


449 




2911 


01:26:04.56 


-37:58:36 


71 


38 


3093 


03:10:55.15 


-47:24:31 


40 




3112 


03:17:58.27 


-44:14:16 


7 


25 


3122 


03:22:18.19 


-41:21:29 


136 




3158 


03:42:54.82 


-53:38:08 


120 


22 


3194 


03:59:08.16 


-30:11:36 


25 




3341 


05:25:32.66 


-31:35:56 


24 




3744 


21:07:13.49 


-25:26:54 


122 




3806 


21:46:27.00 


-57:16:47 


60 


81 


3809 


21:46:57.72 


-43:54:33 


224 




3822 


21:54:09.62 


-57:51:15 


105 




3825 


21:58:27.26 


-60:23:56 


113 




3827 


22:01:55.90 


-59:56:57 


76 





The values of the ellipticities are mostly quite small, 
and even indistinghuisable from 0.0 for quite a few of the 
clusters. This might seem to be in contradiction to other 
results on the elongation of clusters (e.g. Plionis et al. 
1991, and de Theije et al. 1995). However, it must be ap- 
preciated that our ellipticities are apparent ellipticities, i.e. 
not corrected for the effect of the aperture, and that they 
refer to the central regions only (about 5 King core radii: 
i.e. a radius of slightly more than 500 kpc). 

3.2. The values of r c , (3 and ay, 

Using the fitted values of xo, yo, e, and 4> we subsequently 
made MLM fits for r c , (3 and cr?,, for the 62 clusters in Tab. 
1, for each of the 4 profiles, for which the solution of xo, 
yo, e, and <fi converged. 

In Tabs. 4 to 7, we give the values of r c , (3 and for 
all clusters in the sample of 62, i.e. for all clusters with a 
reliable centre position, for each of the model profiles. For 
some of the clusters, one or more of the MLM fits did not 
converge for one or more of the parameters. In those cases 
we give the SIMPLEX values without error estimates. The 
error estimates obtained for each of the 3 parameters, as 
given by MIGRAD are also given in Tabs. 4 to 7, in which 
we also give the values of -21n(L). As is well known, the 
magnitude of this parameter (the maximum likelihood for 
the particular model profile with the best-fit parameters) 



cannot be used to make a statement about the absolute 
probability that a given observation is indeed induced by 
the underlying continuous probability function described 
by the particular profile. 

As a check on the 'meaning' of the error estimates 
from MIGRAD, we have generated about 150000 artifi- 
cial clusters, for which the total number of galaxies, the 
ratio between the number of galaxies in the cluster and 
in the background, and r c - and /3-values for the 4 model 
profiles were chosen in the ranges spanned by the obser- 
vations. To each of the artificial clusters a MLM fit was 
made in the same manner as for the observed clusters. 
The errors in the parameters of the fits to the artifical 
clusters depend somewhat on the type of profile and on 
the assumed parameters, but globally, the results can be 
summarized as follows. 

The background density is recovered with an average 
accuracy of about 25% (with actual errors between about 
10 and 40%), with a slightly better result for the King and 
Hubble profiles than for the NFW and de Vaucouleurs pro- 
files (which is not surprising). Almost identical percent- 
age errors are found for r c . In general, [3 is slightly better 
known, with an average error of about 15% (with actual 
errors between about 5 and 25%). As can be seen from 
Tabs. 4 to 7, the errors given by MIGRAD are essentially 
always in the ranges found for the artificial clusters. 
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The average values of r c and the dispersion around 
the mean arc 128 ± 88, 189 ± 116, 292 ± 191 and 1582 ± 
771 kpc, for the King, Hubble, NFW and de Vaucouleurs 
profiles respectively. There are no general relations for the 
four values of r c that one obtains by fitting an arbitrary 
galaxy distribution with the four model profiles. However, 
from linear regressions between individual values we find, 
in our dataset, that: 
r cK = (0.62 ± 0.07)r cH + (18 ± 53) 
r cK = (0.28 ± 0M)r cNFW + (36 ± 51) 
r cK = (0.06 ± 0.01)r cdeV + (27 ± 70) 

For (3, the average values and dispersions around the 
mean are 1.02 ± 0.08, 1.03 ± 0.07, 0.61 ± 0.05 and 7.6 ± 
0.5, for the King, Hubble, NFW and de Vaucouleurs pro- 
files respectively. Since (3 is closely related to the asymp- 
totic logarithmic slope of the profile, there are predictions 
for the relations between the values of (3 for the four pro- 
files. As can be easily deduced from the four expressions 
in § 3.1, one would expect that (3r = (3k, (3nfw = 0.67 
Pk and PdeV = 8 (3k- It is clear that these relations are, 
to within the errors, obeyed by the observed average val- 
ues for the four profiles. From inspection of Tabs. 4 to 7, 
it is clear that it is not very meaningful to make linear 
regressions between individual /3-values because the dis- 
tributions of (3 for each profile type are quite narrow com- 
pared to the estimated errors in the individual (3- values. 

The values of the average backgrounds are (4.7±2.6), 
(3.7±2.6), (4.0±2.3) and (3.8±1.7) times 10~ 5 galaxies 
per square arcsec for the King, Hubble, NFW and de Vau- 
couleurs profiles respectively. For this parameter, the pre- 
diction is clearly that it should be identical for the four 
fits, as is indeed the case to within the errors. Linear re- 
gressions between the individual values of the background 
density give: 

a bK = (0.89± 0.12)a bH + (-1.6± 1.6)10" 5 . 
a bK = (0.80± 0.ll)a bNFW + (1.5± 1.7)10" 5 
o bK = (0.70± 0.19)<7 6dey + (1.9± 2.3)10" 5 . 

We give also in Tab. 4 the value of the contrast C for 
each cluster. This parameter is defined as the fraction of 
the observed galaxies in a pencil beam that is really in 
the cluster. If N to t is the number of objects of the line of 
sight and if N bc k is the estimated number of background 
galaxies, we have 

C = (N tot - N bck )/N tot 

This parameter has been used in paper IV to select 
the more contrasted clusters, in order to see a narrower 
Fundamental Plane. 

3. 3. The distribution of backgrounds on the sky 

The average backgrounds a b as given in § 3.2 are in very 
good agreement for the four different profiles used, as was 
to be expected. The mean of the four values is about 4 
10~ 5 gals arcsec" 2 , and this mean refers to a magnitude 
limit m bj of 19.5. For the same limit, Colless (1989) pre- 
dicted a background of 3 10 -5 gals arcsec" 2 . Lilly et al. 



(1995) and Crampton et al. (1995) found, from CFRS 
data, 5.6 10 -5 gals arcsec -2 and Arnouts et al. (1996) 
and Bellanger et al. (1995), using Sculptor data, found a 
background density of 2 10~ 5 gals arcsec" 2 for nearly the 
same limiting magnitude. 

These values from the literature vary between about 2 
and 6 10 -5 as do the values of the backgrounds found here. 
Given the large variation in our data for what is supposed 
to be a uniform magnitude limit, one wonders what is the 
cause of this variation, and whether there are correlation 
between backgrounds found in neighbouring directions. 

In Fig. 4, we show the values of a b in galactic coordi- 
nates for the clusters with b<-30°. In the cone thus de- 
fined, we have backgrounds for two-thirds of the ENACS 
clusters with z < 0.1. We have not included the clusters 
with b>-30°, as for those the availability of COSMOS 
data for the ENACS clusters is much less complete. Each 
ENACS pencil beam is indicated by a circle, and the size 
of the circle correlates with the value of a b as found in the 
fit of a King profile. The clusters with b<-30° for which 
no COSMOS data were available, or for which no reli- 
able centre positions could be determined are indicated 
by crosses. 



1=0° . b = -30' 




180° 



Fig. 4. The distribution on the sky of the directions towards 
the 95 clusters with b<-30° (see paper II). For 57 clusters a 
background density is known (circles) and the other 38 clusters 
for which no background density is known are indicated by 
crosses. The size of the circles reflects the background density. 

As can be seen from Fig. 4, there is some structure in 
the distribution of the background values, in the sense that 
large values of a b appear to cluster on scales of about 20 
degrees, which corresponds about to 100 Mpc for z=0.1. 
However, around the south galactic pole, both high and 
low background values occur, and this is also the case 
around 1 ~ 338°, b ~ —46°. The latter concentration of 
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high backgrounds corresponds to a supercluster mentioned 
in paper I; another, at 1 ~ 255°, b ~ —54 corresponds 
to the Horologium- Reticulum supercluster (Lucey et al., 
1983), that was also mentioned in paper I. 

3.4- Which profile fits best? 

The MLM technique does not provide an absolute esti- 
mate of the probability that a distribution of galaxies is 
'produced' by a certain type of profile. However, it does 
allow a comparison of the relative merits of different pro- 
files, through the likelihood ratio statistic (see, e.g. Meyer 
1975). If there are two alternative hypotheses that we are 
testing on a single data-set, the likelihood ratio statistic 
allows one to estimate the probability that the hypothesis 
with the lowest likelihood is false, under the assumption 
that the hypothesis with the highest likelihood is true. 

In practice, one computes the statistic — 2ln(Li/ L2) 
from the values of — 21nL given in Tables 4 to 7, for the 
different density profiles, taking one of the values of the 
likelihood found for a given data-set as a reference L\, 
and comparing that to the likelihood L 2 found for the 
same data-set and a different profile. For large samples, 
the statistic — 21nLi/L 2 has a \ 2 distribution with n de- 
grees of freedom, where n equals the number of parameters 
not restricted by the null hypothesis. 

In the present case, the parameters that are fixed for all 
density profiles (i.e. the centre position, the ellipticity and 
position angle) do not contribute to the degrees of free- 
dom of the distribution, and n — 3 because the core-radius 
(or characteristic radius) r c , the slope (3, and background 
density Oh are estimated independently for the different 
profiles. Note that the central density is not a free param- 
eter, as it is constrained by the normalization implicit in 
the MLM technique. 

In general, the likelihoods obtained from the fitting of 
the King profile are higher than those obtained from the 
fitting of the other profiles: in 37/45, 34/50 and 38/51 
cases, respectively, comparing the King profile with the 
Hubble, NFW and de Vaucouleurs profiles. However, ac- 
cording to the likelihood ratio statistic the differences in 
the likelihood values obtained from fitting different profiles 
to the same cluster data-set, are generally not statistically 
significant (Fig. 5). 

If we set a 95 % significance level for rejection, and take 
the King profile as our null hypothesis, we do not reject the 
Hubble profile as inferior to the King profile for any of the 
clusters. For one cluster, A3528, the King profile provides 
a significantly better fit than the NFW profile, and for 
two clusters (A3128 and A3141) the de Vaucouleurs profile 
is rejected. There is, however, no cluster for which the 
Hubble, the NFW or the de Vaucouleurs profile provides 
a significant better fit than the King profile. 

Most of the individual cluster density profiles, consid- 
ered on the entire available area, do not allow us to select 
one of the four model profiles as the clear favorite and re- 



ject the others. However, in most cases the likelihood of 
the fit is highest when the King model is used. This may 
well be an indication that our inability to choose a partic- 
ular model profile is due more to limited statistics than to 
a fundamental problem with the selection of a best-fitting 
profile. Moreover, we know that the shapes of the four 
profiles are very similar in the outer parts of the clusters, 
beyond 1 or 2 King core radii. The differences between for 
example a King and an NFW profile occur mainly in the 
very central parts of the clusters. 

We have re-computed the statistic — 2\n(Li/ L2) inside 
a square of 2 King core radii for the King and NFW pro- 
files. In general, the differences between the fits are more 
significant in the inner regions: 32 out of 50 clusters show a 
better fit for the King profile at a confidence level of 99%, 
while only 6 out of 50 yield a better fit for the NFW pro- 
file at the same confidence level. This again seems to point 
to the existence of a core in a large fraction of clusters of 
galaxies. 

4. Density profiles of composite clusters 

The result in § 3.4, viz. that the likelihood ratios for the 
King and NFW profile fits, when taken at face value, all 
indicate that the King profile provides a better fit than 
the NFW profile is quite tantalizing, although we realize 
full well that by themselves none of the likelihood ratios 
except one or two really indicate a highly significant pref- 
erence for the King profile over the NFW profile. However, 
it cannot be excluded that the King profile in fact pro- 
vides a better fit for a majority of our clusters, but that 
the limited statistical weight of the data for most of our 
clusters, when taken by themselves, prevents a significant 
demonstration of that fact. The fact that the King profile 
is preferred over the NFW profile in the central regions 
is quite interesting, because this result is obtained with 
smaller numbers of galaxies. 

By combining the galaxy distributions for many clus- 
ters we have increased the statistical weight, in an attempt 
to get a more significant answer to the question if certain 
model profiles fit better than others, and particularly if 
galaxy distributions in clusters have cores or central cusps. 

4-1. The construction of the composite clusters 

The combination of the galaxy distributions of many clus- 
ters to produce what we will refer to as a composite clus- 
ter, requires a lot of care. In order not to smooth (or pro- 
duce) possible cusps, and to avoid artefacts as a result of 
differences in sampling or due to the superposition of elon- 
gated, imperfectly centered galaxy distributions, we have 
used the following procedures. 

First, we must account for the fact that different clus- 
ters have different 'sizes'. If one were certain that all clus- 
ters obey the same type of projected galaxy density profile, 
with different 'core radii' to be sure, one could simply put 
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Fig. 5. Distributions of — 21nLi/L2 for the comparison between King and Hubble profiles (left), between King and NFW profiles 
(centre), and between King and de Vaucouleurs profiles (right). Negative values indicate that the first profile is preferred. We 
also show 95 % (long dash) and 75 % (short dash) significance levels. We note that for the right histogram, A3141 is out of the 
range with a value of -25.72. 



the projected distances between the galaxies in all clus- 
ters on the same scale, by using the core radius to scale 
all projected distances before adding the galaxy distribu- 
tions. However, one can also make a case for scaling by 
r 2 oo, the radius within which the average density is 200 
times the average density in the universe. According to 
Navarro, Frenk and White (1997), scaling with r 2 oo takes 
into account differences in mass. In the present discus- 
sion, we have produced composite clusters without scaling. 
Only for the comparison between the composite clusters, 
i.e. for getting an indication which profile best fits the 
composite clusters (see § 4.4), we have applied scaling. 

A second result of the differences in characteristic scale 
is that the observations of the various clusters do not cover 
the same aperture, when expressed in core radius. To cor- 
rect for this, we have applied a slightly modified version 
of the method devised by Beers and Tonry (1986), and by 
Mcrrifield & Kent (1989). The essence of that method is 
that one adds a model contribution to the observed data 
in the area where a particular cluster does not contain 
data, based on its relative contribution in the area where 
it does contribute. In practice, we combine the galaxy dis- 
tributions of the clusters, beginning with the cluster with 
the largest and second largest core radii. We add some 
galaxies to the second largest cluster in the area where 



this cluster does not have data, according to the density 
profile that was fitted to the central region. For this 'ex- 
trapolation' we used the King-profile fit, but as all model 
profiles that we used have essentially the same outer loga- 
rithmic slope, the result would have been identical had we 
used e.g. the NFW profile fit. This process is repeated for 
all remaining clusters in order of decreasing core radius. 

In this method there is an evident problem with the 
background. When we "reconstruct" the profile as de- 
scribed above, the background galaxies are implicitly, and 
approximately taken into account. The reason for this is 
that the galaxies that are added in the outer regions of 
the smaller clusters, are added in proportion to the total 
number of galaxies, i.e. cluster members plus background 
galaxies. This causes the background to be not exactly 
'conserved'. However, it is not simple or even possible to 
do much better without redshift information. To estimate 
the background effect, we will build a composite cluster 
(see § 4.3) with only ENACS galaxies, for which the mem- 
bership is clear from the redshift. 

4.1.1. The effect of cluster elongation 

We have taken the elongation of the clusters into account 
by 'circularizing' the individual galaxy distributions by in- 
creasing all projected distances orthogonal to the major 
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axis by the appropriate factor, thus 'expanding' the dis- 
tribution parallel to the minor axis. For this correction 
we used the ellipticities from Tab. 1. Even though the el- 
lipticities of most of our clusters are not very large, this 
correction may be important because the superposition of 
elliptical galaxy distributions with randomly distributed 
orientations will cause the outer densities to be underesti- 
mated with respect to the inner ones, which may produce 
an artificial cusp (see Fig. 6). 




Fig. 6. Superposition of 4 elliptical clusters with different ori- 
entations of the major axes. The intensity of the shading is 
proportional to the number of contributing clusters. 

We have checked the importance of the ellipticity cor- 
rection by constructing two artificial composite clusters. 
Both are the superposition of 50 'perfect', artificial King- 
profile clusters, each containing 100 galaxies in a 1200 
kpc square aperture (without background). We assumed 
[3 = 1.0 and selected r c from a uniform distribution in 
the (50,150)-kpc interval. One composite cluster was made 
with an ellipticity distribution which mimics the observed 
one. I.e., we assumed 35 clusters to have e uniformly dis- 
tributed in the interval (0.0,0.1) and the other 15 uni- 
formly in the interval (0.25,0.75). We assumed position 
angles to be randomly distributed between 0° and 360°, 
which is fully consistent with the observations. For the first 
composite cluster, no correction for ellipticity was applied, 
the galaxy distributions were simply summed. 

For the second artificial composite cluster we also used 
50 clusters, again with /3 = 1.0 and r c from a uniform 
distribution in the (50,150)-kpc interval. However, in this 
case, all ellipticities were drawn uniformly from the range 
(0.0,0.1). In Fig. 7, we show the profiles of the two artifi- 
cial composite clusters. As expected, the composite clus- 



ter with the real ellipticity distribution, and without cor- 
rection for ellipticity, has a 2<r-excess in the very centre 
compared to the artificial cluster built from almost round 
clusters. Adding 15 elongated clusters thus induces a small 
but significant cusp. 



-1.2 




1.6 1.8 2 2.2 2. A 

Log t rad 1 us 1 ( kpc ) 



Fig. 7. The two artificial composite clusters, each built from 
50 'perfect' King-profile clusters of 100 galaxies. The profile in- 
dicated by circles is the result of taking the observed ellipticity 
distribution, but without applying a correction for ellipticity. 
The profile indicated by crosses is the result of adding 50 clus- 
ters which are all practically round, making a correction for 
ellipticity unnecessary. Errors are shown as dashed and solid 
lines. 



4.1.2. The effect of centering errors 

Whereas neglecting the ellipticity correction produces an 
artificial cusp, errors in the central position will tend to 
destroy (if not totally, at least partly) a real cusp in a com- 
posite cluster, as argued by Beers and Tonry (1986). The 
magnitude of this effect depends fairly strongly on the ra- 
tio of the position errors and the core radii. We have tested 
the effect of errors in the centre position for our dataset, 
using again 'perfect', artificial NFW-profile clusters. We 
use 50 perfectly circular clusters with 100 galaxies each. 
The characteristic radii are uniformly distributed in the 
range (250,300) kpc. 

We produce five composite clusters, the first with per- 
fect centering and the other with random shifts of the real 
centres lower than 40 kpc, 65 kpc, 85 kpc, 125 kpc and 
with arbitrary orientation. The profiles for the two first ar- 
tificial composite clusters are shown in Fig. 8, from which 
we conclude that there is indeed a smoothing due to the 
position errors, but it is very small. Therefore it is very 
unlikely that our results are seriously influenced by it. Ac- 
tually, a 2-D Kolmogorov-Smirnov test indicates that the 
two artificial composite clusters are indistinguishable at 
a level of about 99%, but this result refers to the whole 
area, and not just the central part of the clusters. In order 
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to quantify the effect of larger shifts, we have compared 
the quality of the NFW and King profile fitting for the 
5 composite clusters through the likelihood ratio statistic 
(see Sec. 3.4). We find that the first 3 clusters (shifts of 0, 
40 and 65 kpc), are significantly better fitted by a NFW 
profile than by a King profile at the 99% confidence level. 
A shift of 85 kpc provides also a better fitting for a NFW 
profile but only at the confidence level of 95%. Finally, for 
shifts of 125 kpc, we have a slightly better fitting for a 
King profile, even if the difference is not significant. This 
shows that even if we attribute the entire difference be- 
tween our isodensity centres and the ROSAT centres to 
errors in our centre positions (which is certainly an over- 
estimation), those errors are not likely to have destroyed 
a potential cusp. 




1.4 1.6 1.8 2 2.2 2.1 
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Fig. 8. Superposition of the profiles calculated for com- 
posite clusters that were built from artificial clusters with 
NFW-profiles with true centres (solid line) , with moderate cen- 
tre shifts (40 kpc: short-dashed line) and large centre shifts 
(125 kpc: long-dashed line). 



4.2. The COSMOS composite cluster 

For the construction of the composite cluster based on 
COSMOS data, we have used the 29 clusters that we also 
used in our analysis of the Fundamental Plane of clusters 
(Adami ct al. 1998). This subset of 29 clusters was selected 
from the sample of 62 clusters in Tab. 1 on the basis of the 
regular character of the galaxy distribution. For these 29 
clusters, we have indicated in Fig. 1 the areas that we used 
in the construction of the composite cluster. It can be seen 
that these are free from evident substructures. In order to 
quantify this fact, we have applied a Dressler-Shectman 
test inside these areas. At a confidence level of 5%, there 
are only 4 clusters among the 29 which show signs of sub- 
structures: viz. A1069, A2644, A3122 and A3128 (14 % of 
all the galaxies). If we exclude the Emission Line Galax- 
ies only A1069, A2644 and A3128 show substructure. The 
contribution of these clusters is certainly minor and we 



considered the composite cluster essentially free from sub- 
structures. A3 128, which is classified as multimodal, is in- 
cluded as its main peak is well-defined. The 29 clusters 
all have a redshift smaller than 0.1 and at least 10 red- 
shifts in the area within 5 King core radii. In principle, the 
latter requirement is not important for the present discus- 
sion, but if one adds the 15 regular clusters from Tab. 1 
with less than 10 redshifts, the total number of galaxies 
increases by only 13%. 

From the 29 clusters with 4735 galaxies, we have con- 
structed a composite cluster; in doing so we added 2505 
galaxies in the outer parts of the smaller clusters (see § 
4.1), to produce a cluster with 7240 COSMOS galaxies. 
The area where we fit the models on the composite clus- 
ter has a radius of 600 kpc, which is similar to that of 
the area over which the profiles of the individual clusters 
were fitted. For the present discussion no scaling of pro- 
jected distances was performed. Although the dispersion 
of the characteristic scales is not very large (see § 3.2), the 
superposition of galaxy distributions with different scales 
may affect the profile of the composite cluster. For ex- 
ample, adding clusters with identical types of profile with 
small and large core radii might produce a cluster which 
does not have the same profile as the individual clusters 
from which is was built. When discussing the question of 
which type of profile best fits the observations (see § 4.4) 
we will therefore also discuss results for composite clus- 
ters for which the constituting clusters were scaled before 
adding them. 

In Fig. 9 we show the result for the COSMOS compos- 
ite cluster, for which all projected distances were scaled 
with the King and NFW characteristic radii of the individ- 
ual clusters: we take the value predicted with the regres- 
sion line between the King and NFW radii. We note here 
that the relation between these two radii is well defined 
(see § 3.2). The profile was calculated in radial bins which 
each contain 20 galaxies. The dashed lines indicate the 1- 
o~ range around the observed values (the dots). The two 
full-drawn curves represent the best King and an NFW 
fits. It is clear that the observed values within r c (log r < 
0) indicate a flat profile. In § 4.4 we will discuss the result 
of a formal test of which of the two model profile best 
represents the three composite clusters. Looking at this 
figure, one may already get some idea about the outcome 
of that test. 

Using the MLM method we fitted both a King and 
an NFW profile to the (unsealed) composite COSMOS 
cluster. These fits give: j3 = 1.00 ± 0.02, r c = 89 ± 5 and 
a b = 2.00 ±0.03 10~ 3 for the King profile, and (5 = 0.56 ± 
0.01, r c = 318±34 and a h = 1.44±0.05 10" 3 for the NFW 
profile. The value of (3 of 0.56±0.01 for the NFW-profile is 
somewhat lower than the value one would predict on the 
basis of the (3 for the King-profile fit, using the expected 
relation between the two (3's. Imposing (3 = 0.67 for the 
NFW profile, the fit is worse but we obtain a b = 1-70 10~ 3 , 
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Fig. 9. The projected density in the composite COSMOS clus- 
ter that was produced by scaling all projected distances with 
the core radii of the individual clusters, so the horizontal scale 
gives r/r c . The dots give the observed values in bins which 
contain 20 galaxies each. The dashed lines indicate the 1-a 
interval on either side of the observed values. The full-drawn 
curves represent the King and an NFW profile. 

which is closer to the value found in the King-profile fit. 
This shows that and at, are correlated. 

The values of r c , and 07, for the composite cluster 
are globally consistent with the average values found for 
the 62 clusters in § 3.2. The agreement is also good if 
we compare with the average values for the subsample of 
the 29 clusters used here rather than for the total sample. 
For the 29 clusters, the average values of and r c and 
the value of £er M are 1.04±0.08, 115±67 kpc, 1.17±0.56 
1CT 3 for the King-profile fits and 0.61±0.05, 276±181 kpc, 
1.13±0.55 1(T 3 for the NFW-profile fits. 

If we restrict the fits to the central region (radius: 300 
kpc) of the composite cluster (924 galaxies), we obtain 
= 1.02 and r c = 106 kpc for the King profile and = 0.55, 
r c = 272 kpc for the NFW profile. These values for the cen- 
tral region agree better with the values for the individual 
clusters; this may be partly due to the fact that the addi- 
tion of galaxies in the outer regions of the clusters hardly 
affect the central region. But the relative importance of 
errors in the assumed centre positions or in the core radii 
determinations is probably larger. 

We have also made a composite cluster from those 
14 clusters for which a central position is available from 
ROSAT, using the latter centre rather than ours. We fit 
both a King and a NFW profile in a more central area 
(radius: 200 kpc) and find a better fit for the King profile 
at the 95% confidence level. The use of the ROSAT cen- 
tres does not change our conclusion about which profile 
fits best. 

4.3. The EN ACS composite cluster 

In principle, the ENACS data could have been used to 
make independent solutions for r c and for the individ- 



ual clusters. Because the ENACS data allow us to select 
cluster members through their rcdshifts, the backgrounds 
would, by definition, be zero. However, the ENACS data 
are limited to the central regions of the clusters, and red- 
shifts are not available down to the magnitude limit of the 
COSMOS data. Therefore, the number of ENACS galaxies 
in a cluster is, in general, too small to make a reliable fit. 
By combining the ENACS data for all 21 clusters among 
the 29 for which the ENACS data cover a rectangular 
area of at least 400 kpc to a side we have constructed an 
ENACS composite cluster with 388 galaxies. 

MLM fits to this composite cluster yield = 1.00±0.05 
and r c = 91 ± 12 for the King profile and = 0.51 ± 0.03, 
r c = 274 ± 61 for the NFW profile. When we impose = 
0.67 for the NFW-profile fit leads to a value of the like- 
lihood that is a factor 5 worse than the maximum likeli- 
hood. It is likely that the lower value of is, at least partly, 
due to the fact that the size of the aperture is fairly small 
compared to the value of r c for the NFW-profile fit. 

Although the statistical weight of the ENACS compos- 
ite cluster is much less than that of the COSMOS com- 
posite cluster, the results of the 2-parameter, rather than 

3- parameter, fit provide strong support for the profile pa- 
rameters found in § 4.2. 

4- 4- Which profile best describes the observations? 

Using the COSMOS composite cluster, we have again ad- 
dressed the question posed in § 3.4. We use three different 
versions of the composite cluster: one without scaling (as 
in § 4.2), one in which all projected distance are scaled 
with the values of r c for the individual clusters as derived 
in § 3.2 (and given in Tabs. 4 to 7), and one in which we 
scaled with the individual values of r 2 oo- The latter were 
calculated from the ENACS velocity dispersions, following 
e.g. Carlberg et al. (1997). For all three composite clus- 
ters we have made fits with a King- and an NFW-profile 
model. 

We thus derived three values for the likelihood ratio 
— 21n(Li/i 2 )- These are -17, -17 and -10 for unsealed, r c - 
scaled and r2oo-scaled composite clusters, respectively. All 
three values are negative, and thus indicate a preference 
for the King profile over the NFW profile. The formal as- 
sociated significance levels are 99, 99 and 95 percent. So, 
at face value, the result of § 3.4 is confirmed and ampli- 
fied: that the majority of our clusters are indeed better 
explained by King profiles than by NFW profiles. 

The following caveats must, however, be remembered 
in connection with this conclusion. 

First: there is a small smoothing effect due to the errors 
in the centre positions, which we cannot quantify very well 
(using the best-fit centres is the best that we can do). Sec- 
ondly, our selection of clusters is certainly not unbiased: 
we have used the most regular clusters (which in general 
are among the most massive ones). This could produce a 
sample of well-evolved clusters, which need not be typical 
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of the total rich cluster population as far as the galaxy 
density profile is concerned. Thirdly, the assumption un- 
derlying the composite clusters, namely that when prop- 
erly scaled , all clusters have the same profile, may well be 
a simplification. Finally, the likelihood ratio for the two 
fits to the composite clusters with r2oo-scaling is least in- 
dicative of a preference for a King profile over an NFW 
profile. This may be just a statistical effect, but it might 
also be an indication that the absence of scaling, or the 
r c -scaling, cause some smoothing due to incorrect com- 
bination of profiles which, intrinsically, have a moderate 
cusp. The likelihood ratio for the r2oo-scaling applies to a 
composite of the 7 clusters in which the data extend out 
to at least 0.75 r 20 o (this cluster contains 1092 COSMOS 
galaxies) . 

4-5. The effect of absolute magnitude 

Finally we calculate likelihood ratios for different ranges of 
absolute magnitude. The reason for this is that e.g. Carl- 
berg et al. (1997) claim the existence of a central cusp for 
the bright galaxies in clusters. We use the COSMOS com- 
posite cluster and we fit King and NFW profiles in the 
central 300 kpc (radius), for different magnitude ranges. 
We define the following 4 intervals of absolute magnitude, 
which contain 200 galaxies each: (-22.4,-18.78), (-18.77,- 
18.03), (-18.02,-17.48) and (-17.47,-16.89). The likelihood 
ratios in these intervals are, +1, -8, -8 and -4, respec- 
tively. It is clear that the preference for a King profile is 
certainly not shared by the brightest galaxies. 

We have tried to estimate the magnitude for which 
the galaxy profile seems to change character. Taking nar- 
rower intervals, which contain 50 galaxies each, we find 
that for absolute magnitudes brighter than -18.4±0.2 the 
likelihood ratio (or rather — 2\n(Li/ L2)) is about 0, while 
fainter than -18.4±0.2 we obtain values between -4 and 
-8. That bright and faint galaxies do not have the same 
distribution is not totally new. E.g., Biviano et al. (1996) 
and Dantas et al. (1997) find that the faint and the bright 
galaxies in the Coma cluster and in A3558 respectively, 
have different distributions. The bright galaxies are very 
concentrated around a few specific positions while the dis- 
tribution of the faint ones is smooth, and without a cusp. 

4-6. Comparison with results in the literature 

From the discussions in § 3.4 and § 4.4 we conclude that 
the COSMOS and ENACS data are better fitted by a King 
profile (or more generally: a profile with a core) than by 
an NFW profile (or rather: a profile with a cusp). In other 
words: our data seem to suggest a logarithmic slope in the 
central cluster region that is flatter than -1. Admittedly, 
this conclusion is based primarily on the composite clus- 
ters, and therefore we cannot be sure that the conclusion is 
true for all of the individual clusters as well. On the other 
hand, the results of the discussion in § 3.4 suggest that 



the conclusion may well be valid for most of the individ- 
ual clusters. In addition, there is some evidence that the 
intrinsically brighter galaxies have a more cusped profile 
than the fainter ones. 

In the literature there are several results that seem 
to be at odds with our first conclusion. Beers and Tonry 
(1986) argue e.g. for the presence of cusps and Carlberg et 
al. (1997) concluded that the CNOC clusters at a redshift 
of ~ 0.3 have galaxy density profiles that are quite con- 
sistent with the NFW profile. Both conclusions are based 
on composite clusters, and as we discussed in § 4.1 several 
details of the construction of a composite cluster influence 
the destruction of cusps or their artificial production. Nei- 
ther Beers and Tonry, nor Carlberg et al. seem to have 
taken the elongation of the clusters into account. As we 
found in § 4.1.1, that may be (partly) responsible for the 
appearance of a small cusp. On the other hand, it is also 
possible (if not likely) that the differences are largely due 
to the fact that our galaxy samples extend to fainter ab- 
solute magnitudes. We note here that the limiting magni- 
tude of the CNOC survey (Carlberg et al., 1997) matches 
well with the value given in § 4.5 for which the King and 
NFW profiles become equivalent. 

Therefore, the apparent disagreement may be largely 
resolved by our second conclusion, viz. that the brighter 
galaxies have a more cusped distribution than the fainter 
ones. Incidentally, this latter conclusion would also seem 
to indicate that for our clusters, the errors in the centre 
positions are indeed sufficiently small that they do not 
erase the cusp in the distribution of the brighter galaxies. 

The values that we find for r c are consistent with 
the results obtained by Girardi et al. (1995) and Bah- 
call (1975). Our value for (3 is only moderately consistent 
with that obtained by Lubin & Bahcall (1994), who found 
0.8±0.1. This discrepancy may be due to the fact that 
these authors fit their profile in the outer regions (500 to 
1500kpc), which may contain substructure. 

Our values of r c and (3 are local (z < 0.1), and our (3 
value can be compared with (3 estimates for other epochs. 
Recently, Lubin and Postman (1996) have studied 79 dis- 
tant clusters from the Palomar Distant Cluster Survey 
(PDCS) with estimated redshifts between 0.2 and 1.2, 
where the precision of the redshift estimation was 0.2. 
They found r c = 50 kpc and (3 = 0.7 using a King profile. 
Because they have estimated the background at 1.2 Mpc 
from the centre, their estimate of r c is probably not very 
accurate, and we do not consider the difference between 
their and our value for r c to be really significant. 

If one accepts the formal error bars for (3, the differ- 
ence in the (3 values for distant and nearby clusters would 
be significant. If the differences between the values of /3 
are real, this indicates a change of the outer slope of the 
density profile with redshift (Tab. 8), in a way that is con- 
sistent with global ideas about an increase of the concen- 
tration of the clusters with cosmological time. According 
to commonly adopted formation and accretion scenarios 
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for clusters of galaxies, young clusters have a flatter den- 
sity profile than older, more evolved clusters. 

5. Profile slope (3 and S7„ 

From the analysis in § 3.2, § 4.2 and § 4.3 it is clear that 
the outer slopes of the galaxy distributions of the individ- 
ual (and composite) clusters are all very much consistent 
with a value of 1.0, for King- and Hubble-profile fits. In ad- 
dition, there is satisfactory consistency with the /3-values 
for NFW- and de Vaucouleurs-profile fits (certainly for the 
individual clusters). We therefore conclude that this result 
is quite robust. 

Quite a few numerical simulations have appeared in 
the literature from which the average present-day slope 
of the density profiles that is predicted as a function of 
the initial fluctuation spectrum, the cosmological param- 
eters etc., can be estimated. In the following, we compare 
our result for (3 with predictions obtained by Walter and 
Klypin (1996), Jing et al. (1995), Crone et al. (1994) and 
by Navarro, Frenk and White (1995). 
Walter and Klypin (1996) simulate a mixed HCDM model 
with VIcdm = 0.6 and f2„ = 0.3. The size of the simula- 
tion box is between 400 and 500 Mpc and they use the 
Harrison-Zeldovit'ch fluctuation spectrum. We will refer 
to it as the HCDMwk model. 

Navarro, Frenk and White (1995, 1996) simulate a CDM 
model with = 1.0. The size of the simulation box is 180 
Mpc. We will refer to it as the CDMnfw model. 
Jing et al. (1995) simulate seven scenarios, in 128 Mpc 
boxes. Each model gives more than 50 clusters which is 
very similar to our real sample. The fluctuation spectrum 
is the Harrison-Zeldovit'ch one. The models are: a stan- 
dard CDM model with n = 1 (SCDM), three flat models 
with O = 0.1, 0.2 and 0.3 (with A = 0.9, 0.8 and 0.7, re- 
spectively) (F101, F102 and F103) and three open models 
with Q = 0.1, 0.2 and 0.3, and A = (OpOl, Op02 and 
Op03). 

Crone et al. (1994) simulate three models with f2 = 0.1 
and 0.2 using different fluctuation spectra: P(k)ock°, k and 
k 2 . The models are three standard CDM ones (SCDMcO, 
SCDMcl and SCDMc2), three flat universes with fl = 
0.2 and A = 0.8 (F102c0, F102cl and F102c2), three open 
universes with Q = 0.1 (OpOlcO, OpOlcl and Op01c2) 
and three open universes with Q = 0.2 (Op02c0, Op02cl 
and Op02c2). 

All four groups fit 3-D power laws r~ 2/33D to the outer 
regions of the clusters, at the present-epoch which cor- 
responds very well to the redshift range of our clusters. 
Their values are transformed to 2-D values, as follows: 
/?2D = 2133 2 1 ■ W e compare these 2-D model values with 
our observed mean /3 2 d = 1-02 ± 0.08 for the King-profile 
fit in Fig. 10. 

We find that only fairly low values of flo produce aver- 
age values for f3 2 D that are similar to the observed value; 
this is in agreement with other recent results on Q - How- 



ever, not all low-fio models produce the observed value of 
@2D- Taking the results in Fig. 10 at face value, it would 
seem that not only needs fio to be small (a few tenths), 
but it is also not likely that we need to invoke a flat Uni- 
verse (i.e. a non-zero value for A). The only exception to 
this general conclusion seems to be the flat model with Q 
= 0.2 and a fluctuation spectrum P(k)ock°. 

It should be realized that the conclusion that /?2_d is 
best predicted by low-^o models (and possibly mostly 
open models ) is based on a rather incomplete range of 
models. In other words: we cannot exclude models that 
have not yet been calculated or published. E.g. it is too 
early to firmly exclude HCDM models with other fluc- 
tuation spectra than the one calculated by Walter and 
Klypin. On the other hand, the data in Fig. 10 seem to 
suggest that for (3 the influence of the fluctuation spec- 
trum is quite limited, so our conclusion of a low-f^o model 
(most likely open) is probably quite robust. 

This value of flo from the outer slope of the projected 
galaxy distribution can be compared to direct estimates of 
f2 based on cluster masses and luminosities and the field 
luminosity density. The idea is quite simple: one estimates 
the cosmological volume which contains the same mass as 
that contained in the cluster, from a comparison of the 
luminosity density in the field and the total luminosity of 
the galaxies in the cluster. This yields the average mass 
density which yields fig through division by the critical 
density p c . In practice, the method amounts to calculating 
the ratio of the ratio in cluster to the critical ^ ratio in 
the field, and this result does not depend on H a . A recent 
application of the method was discussed by Carlberg et 
al. (1996) for their sample of 16 clusters at a redshift of « 
0.3. 

We have determined luminosities and velocity disper- 
sions for our sample of 29 ENACS clusters discussed in 
paper IV. The luminosities and dispersions are calculated 
in an aperture of five core radii. We derived the projected 
virial masses (e.g. Perea et al. 1990) and obtained the ^ 
ratios for the 29 clusters. The values of of the 29 clus- 

ters range from about 100 to 1000; the average value is 
454, the median is 390 and the robust bi-weight estimate 
of the mean is 420, in good agreement with the value given 
by e.g. Bahcall et al. (1995). The local critical ^£ ratio in 
the field was determined by e.g. Efstathiou et al. (1988) 
and Loveday et al. (1992), who find a best estimate of 1500 
w ith the value most probably in the range 1100 to 
2200. This yields an estimate of fio on the basis of our 
sample of 29 local ENACS clusters of 0.28 ± 0.19. This 
value is thus quite consistent with the low Q value that 
we obtained from the outer slopes of the projected galaxy 
distributions. 

The uncertainty in the 'dynamical' f2o estimate is to 
a large extent due to the large spread in ^ ratios for 
our clusters. It appears that the individual values of ^£ 
correlate moderately well with velocity dispersion, in ac- 
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cordance with the result in paper IV. Expressed in fl we 
find: n = (3.9±1.1) 1CT 2 a v /l00 + (8±157)1(T 3 . If we 
use the relation between luminosity, scale factor and veloc- 
ity dispersion that we derived in paper IV we obtain: Qq 
= (3.2±0.5) 1CT 4 a x v m R-° 19 , which is totally consistent. 
It thus appears that there is a significant dependence on 
cluster velocity dispersion in the determination of SIq from 
cluster ™ ratios, which could easily produce a serious bias 
towards high values. 

Of course, we have assumed in this analysis that in 
clusters light traces mass. Actually, this assumption is not 
really demonstrated, but is commonly used, for example 
in Carlberg et al. (1996). One might wonder how this as- 
sumption could affect our last firj estimate. However, for 
the low a v clusters, the present estimate is in good agree- 
ment with that based on the slope of the density profiles, 
and that gives some confidence that a possible bias is prob- 
ably quite small. 

6. Conclusions 

We have studied the projected galaxy distributions in 77 
clusters from the ESO Nearby Abell Cluster Survey. The 
present sample is an unbiased subset of the volume- limited 
ENACS sample, and thus forms a representative local (z 



< 0.1) sample of rich (Raco > 1), optically selected clus- 
ters. We used both COSMOS and ENACS data to test the 
character of the projected galaxy distributions. In particu- 
lar, we have investigated whether the galaxy distributions 
in rich clusters have cusps or cores in their central regions. 

We have made maximum Likelihood fits to the ob- 
served distribution of COSMOS galaxies to solve for the 
position and the elongation of the clusters. For 15 of the 77 
clusters, no reliable centre could be determined and these 
clusters were not considered further. Using the positions 
and elongations, we subsequently solved for each of the 
62 remaining clusters the three parameters that describe 
each of the four theoretical profiles that we tested, as well 
as the density of background galaxies. The four model pro- 
files that we tested against the data are the King, Hubble, 
NFW (Navarro, Frenk and White) and the de Vaucouleurs 
profiles. Although the solutions do not converge for all of 
the clusters nor for all four profiles, we obtain reliable re- 
sults for between 75 and 95 % of the clusters (depending 
on the model profile). 

We find mean values for r c , the characteristic scale of 
the 2-D galaxy distribution, and dispersions around the 
means of 128 ± 88, 189 ± 116, 292 ± 191 and 1582 ± 
771 kpc, for the King, Hubble, NFW and de Vaucouleurs 
profiles respectively. The outer logarithmic slopes of the 
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distributions were generalized by the usual /3-paramctcr, 
which we find to have the following average values: 1.02 ± 
0.08, 1.03 ± 0.07, 0.61 ± 0.05 and 7.6 ± 0.5, for the King, 
Hubble, NFW and de Vaucouleurs profiles respectively, 
which are consistent. The average background density at 
the limit of the COSMOS catalogue is about 4 10~ 5 galax- 
ies arcsec -2 . 

In order to investigate whether the galaxy distribu- 
tions in our clusters preponderantly have cores or cusps, 
we have determined the likelihood ratio for the King and 
NFW profiles. Using all galaxies down to the COSMOS 
magnitude limit of about m;,. w 19.5, we find that in gen- 
eral the King profile is more likely to be a good repre- 
sentation of the data than the NFW profile. However, for 
the individual clusters this preference for the King profile 
is generally not statistically significant. If we restrict the 
analysis to the central regions, the significance of the pref- 
erence for the King-profile fits increases, even though the 
number of galaxies decreases. 

We have increased the statistical weight for the likeli- 
hood analysis by combining the galaxy distributions in a 
subset of 29 of the 62 clusters, which show a regular galaxy 
distribution. We take special care to avoid the creation of 
an artificial cusp (by taking the ellipticities into account), 
and to avoid the destruction of a real cusp by summing 
distributions with different scale lengths. We have also 
checked that it is unlikely that the uncertainty in the cen- 
tre positions has erased a cusp. For the test we summed 
without scaling projected distances, after scaling with r c , 
as well as with r2oo- In all three cases we find that the King 
profile provides a better fit to the data than the NFW pro- 
file, at confidence levels of more than 95 %. Interestingly, 
this preference is not shared by the brighter galaxies. 

Finally, we have used the outer profile slope (i.e. the 
result that (3 is very close to 1.0), in combination with sev- 
eral results from numerical models to conclude that the 
density parameter Qq is likely to be considerably smaller 
than unity. In addition, the available models indicate that 
the Universe probably has an open geometry (i.e. no clo- 
sure through A is indicated). This low implied value of Qq 
is fully consistent with a direct determination based on 
the ratios of our clusters. 
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VII: Density profiles of rich clusters of galaxies 



Table 4. Fitted parameters for the King models and contrast 
(last column). 
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Table 5. Fitted parameters for the Hubble models. 
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Table 6. Fitted parameters for the NFW models. 
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Table 7. Fitted parameters for the de Vaucouleurs models. 
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Table 8. (3 values for different redshifts from the present study 
and the study of Lubin & Postman (1996). 



redshift 0.07±0.05 0.3±0.1 0.6±0.1 1±0.2 
P 1.02±0.08 0.73±0.17 0.74±0.17 0.68± 0.40 



